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Abstract. Recently, it was shown that an infinite perfectly conducting (PEC) cylinder can 
be nearly perfectly cloaked from normally incident electromagnetic waves using a single- 
layer homogeneous dielectric cladding. Here we study the electromagnetic response of 
such structures with the goal to understand the main mechanisms underpinning this cloaking 
phenomenon. We introduce a simple model of the cloaked PEC cylinder, replacing it by 
an omnidirectional electric-line scatterer and a bipolar magnetic one; accordingly, the far 
field is determined in a compact closed form. The analysis of the results shows that the 
optimal cloaking regime corresponds to the frequency point where the total electric moment is 
drastically mitigated and thus the radiation pattern of the device resembles that of a magnetic - 
dipole line. In the vicinity of the optimal cloaking frequency we observe the response close to 
that of Huygens' pairs of electric and magnetic line scatterers. 
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1. Introduction 

Cloaking devices render objects partially or wholly invisible to parts of the electromagnetic 
spectrum by drastically reducing the total scattering cross section of the objects. Due to 
the obvious interest in potential applications and due to their intriguing nature, cloaks have 
attracted considerable attention from both practical experts and theoretical scientists during 
the previous decade. In the pioneering works HH21, the control of electromagnetic fields 
in order to not to interact with an object is theoretically formulated. They use optical 
conformal mapping techniques and complex coordinate transformations which require exotic 
inhomogeneous and anisotropic materials. On the other hand, actual structures realizing 
broadband electromagnetic cloaking with, e.g., transmission lines, are presented in the papers 
E HI, where experimental validation is provided. Furthermore, the so-called "scattering 
cancellation" technique has been introduced in 0, where drastic reduction of the scattering 
response of a specific dielectric object is accomplished with a dielectric covering layer. 

In this paper we consider cloaking of infinite perfectly-conducting (PEC) cylinders using 
the simplest possible configuration: a uniform dielectric cover O. A PEC circular cylinder 
is covered by a single cladding and is excited by a plane wave (whose electric field is parallel 
to the cylinder axis), formulating a two-dimensional (2D) boundary-value problem. This 
geometric configuration is the same as that for scattering cancellation cloaking of dielectric 
cylinders [5], but the physical phenomena leading the the cloaking effect are different and the 
known design methods for scattering-cancellation cloaks cannot be used for this polarization 
0. In fact, the equations for the required thickness and material parameters of the cloaking 
cladding have no solution for the case of a PEC cylinder in the center [|5l|7]|. However, similar 
devices have been analyzed in electrostatic [8], in quasi-static [0 and in electrodynamic (61 
regimes and satisfying cloaking performances have been achieved. It is therefore important 
to understand the phenomenon and develop a predictive analytical model, enabling the design 
of cloaks for elongated conductive objects. 

Here we introduce an approximate model that can be used in explaining and interpreting 
the aforementioned interesting behavior. In particular, we note that the scattering pattern of 
cloaked PEC cylinders in the vicinity of the optimal cloaking frequency is similar to simple 
dipolar objects. In addition, we expect that the optimal cloaking regime should correspond 
either to minimization of induced lowest-order moment(s) or to the electric and magnetic 
moments forming a Huygens' pair with the minimum scattering in the forward direction. This 
brings us to the idea to model the cylindrical volume by one electric and one magnetic dipole 
line located along the central axis. The corresponding moments are evaluated analytically 
from the canonical field solution and the scattered field of the equivalent structure is expressed 
in a compact form. After validating the model by comparing the data originating from the 
exact solution with that given by the approximate model, we discuss the frequency response 
of the equivalent system. Certain conclusions related to the principle of the device operation 
and justification of the cloaking phenomenon are drawn. Some preliminary results of this 
study have been reported in a conference fTTOll . 
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2. Electromagnetic Fields 

Let us consider a simple two-dimensional (2D) structure depicted in Fig. [TJ An infinite 
PEC cylinder of the radius g is covered by a single isotropic and uniform dielectric layer 
of the thickness (a — g). The cladding layer (Region 1) is filled with a magnetically inert 
material of the relative permittivity e r . The whole structure is placed in vacuum (Region 0, 
the permittivity e , and the permeability /i ). The axis of the cylinder coincides with the z 
axis of the cylindrical coordinate system (p, (p, z), which can be used interchangeably with 
the corresponding Cartesian one (x,y,z). This cylindrical configuration is illuminated by 
a z-polarized (TE or TM 2 in an alternative definition) plane wave of unitary electric-field 
magnitude propagating along the positive x semi-axis ffTTTl : 

+ 0O 

E 0iinc = zE z0iinc (p, (j)) = z ^ A {n)J n {k o p) cos(n0), (1) 

n=0 

where A (n) = j~ n , k = 27r/ x /e /xo is the free-space wavenumber and / is the 
operational frequency. The notation J n concerns the Bessel function of n-th order, while 
the symbol S n0 corresponds to the Kronecker delta. The harmonic time dependence is of the 
form e +j27r ^' and suppressed throughout the analysis. 

Due to the 2D nature of the studied configuration and the isotropy of the materials, the 
(secondary) electric field in each of the two regions (0, 1) is z-polarized and can be expressed 
in the following forms: 

+0O 

E z0 ,scat{p, <P) = J2 B o(n)Hi 2) (k p) cos(n0), (2) 

n=0 

E zl (p,<f>) = [A 1 {n)J n {kp)+B 1 {n)H^{kp)] cos(n0), (3) 

n=0 

where Hn^ = J n — jY n is the Hankel function of n-th order. The wavenumber in region is 
defined as k = k^^fe r . 

By imposing the necessary boundary conditions around the cylindrical surfaces p = g,a, 
the unknown coefficient sequences B (n), Ai(n), Bi(n) of §2§, ^ are readily determined. 
The azimuthal magnetic component in region 1 is expressed as: 

H*iM) = -T^kY" [A^J'Jkp) + B x {n)H^'{kp)\ cos(n0), (4) 



where Co = yPo/ e o — 1207T VI is the wave impedance of vacuum. Based on these field 
quantities, the equivalent dipole moments can be rigorously evaluated. Furthermore, the far- 
field response of the cylindrical structure is determined with the use of Hankel's expansion 
for large arguments: Hn\k p) ~ J ^—e~^ kop ^~^ :[ \ k Q p — > +oo. In particular, 



2 

E z0 ,scat(pA) ~ J^e-^-5>^5 o (n)j n cos(n0), hp +oo. (5) 
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3. Electromagnetic Moments 

3.1. Equivalent Current 

The surface current induced at the PEC surface p = g is given by: 

K(0) = zF 01 (<7,0). (6) 

By definition, the polarization current over the cross sections of the cladding layer is written 
as follows: 

k 

J P oi{p, 0) = zj-^(e r - l)E z i(p, 0). (7) 

Consequently, the total induced current which is the source of the scattering field by the 
cylindrical structure possesses the form: 

J(p, 0) = K(0)5(p -g) + J pol {p, 0) 

= zH^(g, 4>)6(p -g) + z Ae r - l)E zl (p, 0), (8) 

so 

where 5 is the Dirac function. 



3.2. Electric Moment 

The electric moment p of electric current density J existing in a volume (V) is defined in lfT2l 
and can be written in the 2D case as an integral over the cross-section surface (S): 

1 1,1 J dV p* = — (( JdS, (9) 



j^fJJJ(v) 3 2 ^fjJ(S) 
where p* is the electric moment per unit length (p.u.l.) of the z axis. If one substitutes 
expression ([8]) into formula ([9]), only the 0-th term from the harmonic series Q, Q survives; 
accordingly, the following expression is obtained: 

2tt f kl{e r -l){AMVAg,a,k) + BMV H {g,a,k)} 



z- 



k 2 ( c I -kg A 1 (0)J{ ) (kg) + B^H™ {kg) 



(10) 



where c = -J= is the speed of light. The integrals of the Bessel and Hankel functions are 



analytically evaluated: 



Vj{x,^,k)=l J {kp)pdp= , (11) 



x 



V H (x,y,k) = J H o {kpjpdp = . (12) 

This quantity p* = z p* z (in Coulomb) expresses mainly how powerful is the total induced 
electric current for given device dimensions and the frequency. It can be considered as a 2D 
electric dipole line (or an electric current line) along the central axis of the structure, which 
models the electric response of the whole cylindrical configuration. 
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3.3. Magnetic Moment 



The magnetic moment m of a volumetric distribution of electric current density J existing in 
volume (V) is defined in lfT3l and can be rewritten in the 2D case as a surface integral over 
the cross section (S): 

m= lJJJ { rx ^ d V^ m * = lJ^ rxJdS, (13) 

where r = pp + zz is the observation vector in cylindrical coordinates and m* is the magnetic 
moment per unit length (p.u.l.) of the z axis. Obviously, the vector (r x J) has exclusively 
an azimuthal (<fi) component. Due to the orthogonality of the cylindrical functions, only the 
1-st term from the harmonic series ([3]), ([4]) survives. Writing the unit vector along in the 
Cartesian coordinates: <f> = — x sin + y cos 0, we see that the vector of the magnetic dipole 
moment is parallel to the y axis (the term x sin gives an odd function with zero contribution 
to the integral). Accordingly, one obtains the result: 

_ jir f k 2 (e r - l)iA 1 (l)Wj(g,a,k) + B 1 (l)W H (g,a,k)} 

" y 2^Co I -kg 2 k(l) J[(kg) + B 1 (l)H[ 2) '(kg) 



m 



.(14) 



Again, the related integrals of the Bessel and Hankel functions are rigorously determined as 
follows: 



J\[kp)p 2 dp 



jj 2 J 2 (kij)-x 2 J2(kx) 
k 



(15) 



W H (X,M= / H[ 2 \kp)p 2 dp = Wj(x^,k) 



J_ 
2 



^ 3 M 



- X 3 M 





1 




1 


1 


( 


2 




2 


2 




-1 




3 
2 


-1 





1 




1 


1 


( 


2 




2 


2 




-1 


• 


3 
2 


-1 



kip 1 
' 2 ' 2 



2 ' 2 



(16) 



where M is the Meijer function iTPfll . This quantity m* = y m* (in Ampere-meter) expresses 
mainly how powerful is the spatial variation of the induced electric current over the cross 
section of the device for given device dimensions and the frequency. It is a 2D magnetic 
dipole line positioned along the central axis of the structure. 

In the following we will see that this simple model adequately represents the response of 
the cloaked cylinder near the optimal operational frequency despite the fact that the variations 
of the fields inside the cloak are quite significant (its electrical size is considerable due to high 
permittivity of the cladding cylinder). 

3.4. Radiation Fields 



The radiated electric field E p by an electric moment p* p.u.l. of the axis z in a 2D configuration 
is written as follows [fT5l : 

1 



E p (p,0) = -Vx [VG(p, 



x P 



(17) 
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where G is the scalar Green's function in two dimensions, expressed in cylindrical 
coordinates. It corresponds to the field at a point (p, 0) developed from a singular source 
at (P, $). In our case, the position of this dipole coincides with the z axis (P = 0) of the 
coordinate system; therefore, the respective expression takes the form lfT6ll : 



G(p,0) = -^ 2) (M- (18) 

Similarly, the radiated electric field E m produced by a magnetic-moment line m* p.u.l. at the 
axis z in a 2D configuration is given by [fT5l : 

E m (p, 0) = -j27r/p VG(p, 0) x m*. (19) 

If one takes into account the polarization of the moments (p* = zp*,m* = y m*), 



substitutes Green's function ( [18] ) and adds the two electric fields ( [17] ), ( fT9[ ), the total developed 
field (E p + E m ) from the moments couple possesses a single z component given by: 



Ko^at(pA) = % \ m * v H[ 2 \k p) cos — jcp* z HQ 2 \k p) 



4 

The corresponding far- field expression reads 



(20) 



The two electric fields E z0!Scat (p, 0) and E' z0scat (p, 0) represent the same physical quantity, 
namely the field scattered by the cylindrical structure under the plane wave excitation. The 
difference is that the first one (the unprimed, denoted by: E z0jScat (p, 0)) is evaluated from the 
explicit forms of the fields determined through the rigorous canonical solution of the scattering 
problem, while the second one (the primed, denoted by: E' z0 scat (p, 0)) is computed from the 
approximate model of two scattering dipole lines. 

4. Total Scattering Quantities 

By repeating the same procedures described above, one can readily find the corresponding 
field quantities E z0tScat (p, 0), E' z0 scat (p, 0) which describe the situation of a bare PEC cylinder 
of radius g in the absence of the cladding (our solution for e r = 1). In this way, one can define 
the normalized total scattering widths in each case as follows: 

fo n \E zO)Scat (p,(p)\ 2 d0 

CTnorm = llHl \ ~ ; — , (22) 

y \ E Wscat{p^)\ 2 d(t> .... 

a norm = llm ^ ~ ' V 23 ) 

lo \Ko,scat(P,<f>W 

These quantities measure the reduction of total scattering from the PEC cylinder due to 
cloaking cladding (when a norm , cr' norm — >■ 0, we have perfect cloaking effect). Also here 
Cnorm and (j' norm give the normalized scattering widths computed from the rigorous solution 
and the dipole moments, respectively. 

The forward scattering theorem [fT71 is a well-known lemma, frequently used in electro- 
magnetics, which relates the scattered field by a lossless object E z0 scat (p, 0), E' z0 scat (p, 0), 
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along the forward direction (namely the one of the incident plane wave, here: = 0) with the 
total scattered power from the object. In particular, this power (expressed in Watt) is given 
by: 

V2 



scat 



PL 



scat 



hCo 
2V2 



+00 



n=0 



B (n)f 
[cp z - m y \ , 



(24) 



(25) 



respectively to which method we adopt (the full- wave solution or the moments model). In 
other words, in the case of a perfect cloak, the electric quantity £y[cp*] should be equal and 
opposite to the magnetic quantity m*]. The two sources ^[cp*], $s[— m*] are cooperating 
with each other to produce mutually neutralized scattering. 



5. Numerical Results 



5.1. Parameters 

Before proceeding with the discussion of the results, we need to define the value ranges of 
the problem parameters. As the reference operational frequency of the cloak, we assume 
/o = 3 • 10 8 Hz corresponding to unitary free-space wavelength A = 1 m. The physical size 
of the considered structure is kept constant throughout the numerical simulations to simplify 
the procedure and highlight the effect of variations of the other parameters. More specifically, 
we take g = 0.05A and a = 0.08A . This assumption does not affect the generality of 
our investigation because similar dependencies are observed for other dimensions. Since 
the operational frequency does not differ significantly from f , the electrical dimension of 
the PEC cylinder is close to 10% of the free-space operational wavelength. In other words, 
cloaking (even not perfectly) this conducting volume is not a trivial task since its electrical 
size is far from infinitesimal even in terms of the free-space wavelength. 



5.2. Graphs 

In Fig.|2(a)|we sweep the relative permittivity of the cladding e r and evaluate the normalized 



total scattering width using the exact field formulas. The operating frequency is kept constant 
f = fo and the optimal cloaking result (the minimum a norm ) is obtained for e r = 60. 
Therefore, the dielectric constant is kept fixed to this value hereinafter to emphasize the 
frequency response of the device. Needless to say that additional minima are recorded 
for higher frequencies, but we are interested in the cloaking regime at the lowest possible 
frequency (not so high permittivity and electrical size). In Fig. |2(b)| we plot the same 
quantity a norm with respect to the frequency. The optimal frequency f opt is very close to 
fo (fopt = 0.992/q) but not exactly equal to that since the choice e r = 60 does not correspond 
to the optimal case of Fig. 2(a)| Additionally, the shape of the two curves in Figs.[2]are similar 



since the horizontal axis corresponds to a common quantity: the effective electrical distance 
inside the cladding. 
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In Figs. [3j we "zoom" in the operation of the device close to the optimal frequency 
/ = fo-pt to better understand the cloaking effect. In particular, we assemble five graphs of the 
normalized scattering patterns of the cylindrical structure, defined by \imk p^+oo |^ z ° ,scat ^'^| 
for f/f pt = 0.95,0.98, 1, 1.02, 1.05. We observe that the polar plots differ substantially 
from each other despite the smallness of variations of /. More specifically, a backward 
scattering for / < f opt and a forward scattering effect for / > f opt has been recorded. 
Exactly at / = f opt , we obtain a bipolar plot which means that the first (n = 1) term in 
series Q, ([5]) becomes dominant in defining the scattering pattern. To put it alternatively, the 
omnidirectional zeroth term (n = 0) is suppressed at / = f opt and increases rapidly when 
going away from this point. Such behavior reminds us clearly of the scattering cancellation 
principle [5] which is reported as being not feasible for PEC cores and the considered TE 
(referred to as TM 2 in flVj) polarization. Consequently, the present cloaking technique can be 
characterized as an "expanded scattering cancellation" suitable for cloaking of PEC cylinders. 
In fact, we have achieved a de facto scattering cancellation effect for PEC objects. Such a 
feature has been briefly referred in |fl~9ll . 

To study the phenomena in the vicinity of the optimal cloaking frequency / = f opt , we 
adopt the simplified analytical approach with the equivalent moments described in Section III. 
We represent the variation of the normalized total scattering width when computed with the 
exact field formulas (a norrn ) and when evaluated from the respective moments (cr' norm ) in 
the same diagram (Fig. |4j), as functions of the normalized frequency f / fopt- The former 
graph is the same as Fig. 2(b)[ and the latter exhibits remarkable agreement with the first one, 
especially for smaller frequencies. In the vicinity of the cloaking regime, the model certainly 
"captures" the physics in terms of the minimized scattering; implicitly, optimal frequency 
f' t in the moments case is only slightly lower than f opt (f t = 0.984/ ). However, the 
exact level of the maximum reduction cannot be predicted successfully since o' norm is much 
smaller than the actual o norm . This difference is due to the fact that the higher-order terms 
in > 1) in d2|), (Bl) are significant at / = f opt , f' opt where the scattering due to the leading 



term is mitigated. Their effect is mutually canceled in the radiation pattern of Fig. 3(c) but 
they indeed substantially contribute to the overall scattered power and define this gap between 
the levels of a norm and cr' norm as shown in Fig.[4j Moreover, when the operational frequency 
gets higher, the moments model becomes inaccurate since it replaces the entire structure by 
a pair of dipoles. In other words, electrically thick cylinders, that imply strong azimuthal 
variations of induced polarization, cannot be fitted with the few degrees of freedom of our 
approximate moments model. However, the model results exhibit satisfying similarity with 
the actual field waveforms at the moderate frequencies and imitate successfully the developed 
physical mechanism at the cloaking regime. 

Let us now analyze the Figs. [5j which are analogous to Figs. |3^>ut with the use of the 
moments model (the quantity represented is defined as: limfe 0P ^ +oo ,sca ^ P '^| ). Again, the 
cloaking effect is observed during the transition from the backward to forward scattering 
in the radiation pattern of the structure. Exactly on the optimal frequency f t the coated 
cylinder radiates as a magnetic dipole line, namely only the term proportional to m* in (21 ) 
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is of significant magnitude. In this way, we reach the conclusion that the other term of (21 ) 
corresponding to the electric moment cp* z is suppressed. In this way, it becomes evident that 
our model combines the electric and magnetic effect into one formula and assigns to each of 
them a different azimuthal profile: omnidirectional and bipolar respectively. As reported in 
Fig. [4] the values of the scattering response are much lower than in the exact case. 

After validation of the model in the vicinity of the cloaking frequency f t = f opt , 
we can make exclusive use of it hereinafter. In Fig. [6] we show the magnitude variation of 
the two amplitudes incorporated in (21 ). The first one (c|p*|), corresponding to the electric 
response of the structure, is minimized (almost nullified) at / = f t as predicted above. The 
second one (|m*|), determining the magnetic response of the cylinder, increases smoothly 
in the neighborhood of f t . Again, it is clearly noted that the electric moment represents 
the omnidirectional term of the far-field pattern, while the magnetic moment is proportional 
to the bipolar term as defined by expression (21 ). Therefore, the findings from the graphs 
are in accordance with those of Figs. [3] (based on the explicit field formulas), since the 
omnidirectional component is mitigated at the cloaking frequency. 

In Figs.|7J we show the variations of the constituent parts of the quantities cp* z and (— m*) 
(real in Fig. 7(a)| and imaginary in Fig. |7(b)[) with respect to the normalized operational 



frequency f/fLf The negative sign in the magnetic moment is compatible with (25 ), where 
the scattering power is written as the result of the action of two sources: one with magnitude 
Q^ep*] and the other with magnitude 9f[— m*). From Fig. 7(b) it is apparent that the two 
quantities cannot give a zero sum at any frequency, since the imaginary parts (responsible for 
scattering losses) retain the same sign (negative) at all frequencies. Therefore, it is sensible for 
the optimal choice to be / = f t since: (i) both imaginary parts are small as demonstrated 
in Fig. 7(b) (ii) the real part 9ft [cp*] is nullified at / = f t and its magnitude increases 
rapidly slightly far from it as shown in Fig. 7(a) and (hi) the real part 3ft[— m*) does not vary 
considerably in the vicinity of / = f' opt as again indicated in Fig. 7(a) When observing 
the frequency dispersion of the real and the imaginary parts of the electric moment, one 
concludes that it exhibits a typical plasma dependence tending to large negative values at 
low frequencies. On the other hand, the dispersion of the magnetic moment resembles a 
Lorentz-type curve with a resonance at higher frequencies. 

In Figs. [8} we represent the real and the imaginary parts of the scattering field in the 
forward direction (0 = 0) in the far region as functions of / / f' opt from expressions (|5j), (21 ), 

where the common factor y^j^ > e ~^ koP ~^ has been dropped. By taking into account that 
the scattered power by the entire cylindrical structure is proportional to the real part of the 
electric far field as stated in (25 ) and observing Fig. 8(a)[ one concludes again that the perfect 
cloaking is not possible and only a minimization of the overall response is feasible. As far 
as the imaginary part of the forward-scattering quantity is concerned, it is kept moderate for 
/ = f' opt regardless of the method (exact or approximate) one employs to evaluate it. 



On the Minimal Scattering Response of PEC Cylinders in a Dielectric Cloak 



10 



6. Conclusions 

An analytical model based on the equivalent moments of a cylindrical PEC rod covered by 
a conventional dielectric is used to analyze the minimal scattering response of the structure 
when it is illuminated by a normally incident plane wave (with electric field parallel to the 
cylinder's axis). The adopted approximate model is reliable for moderate electrical sizes 
and captures the scattering effect of the illuminated object. In the vicinity of the cloaking 
frequency, a drastic change in the radiation pattern of the covered PEC cylinder is observed; 
in particular, there is a clear transition from backward to forward scattering. It has been found 
that the minimum scattering is achieved when the electric moment is mitigated and in that 
case the system radiates as a magnetic dipole line. 
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Figure 1. A simple 2D structure of a PEC cylinder of radius g covered by a magnetically inert 
layer of thickness (a — g) and relative permittivity e r , is illuminated by a plane wave. 




Figure 2. The normalized total scattering width o norm as function of: (a) the cladding relative 
permittivity e r (/ = /o) and (b) the normalized operating frequency f / f op t i^r — 60). Plot 
parameters: / = c/A = 3 • 10 8 Hz, f opt = 0.992/n, g/X = 0.05, a/X = 0.08. 
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Figure 3. The normalized radiation patterns lim fcri0 _j. +oo !f *°' 3Ct "^ p 'tl! of the cylindrical 
structure determined from the exact field solution for: (a) / = 0.95/ opt , (b) / = 0.98/ opt , (c) 
/ = fopt, (d) / = 1.02/^ and (e) / = 1.05/ opt . Plot parameters: /„ = c/A = 3 • 10 s Hz, 
f opt S* 0.992/q, g/\ = 0.05, a/X = 0.08, e r = 60. 
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Figure 4. The normalized total scattering width as function of the normalized operating 
frequency / / f opt when evaluated through the exact field solution (u norm ) and when computed 
from the approximate moments model (cr' norm ). Plot parameters: / = c/Xq = 3 • 10 8 Hz, 
f opt = 0.992/q, g/X = 0.05, a/X = 0.08, e r = 60. 
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(d) (e) 
Figure 5. The normalized radiation patterns limfe QP ^ +00 |^ 0,3cat | p '^|| of the cylindrical 

I zO, scat vr' V J I 

structure determined from the approximate moments model for: (a) / = 0.95f' opt , (b) 
/ = 0.98/^, (c) / = f opt , (d) / = 1.02/^ and (e) / = 1.05/^. Plot parameters: 
fo = c/A = 3 • 10 s Hz, f t = 0.984/q, g/X = 0.05, a/X = 0.08, e r = 60. 
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Figure 6. The magnitudes of the moments per unit length (p.u.l.) of axis z (c p*,m*) as 
functions of the normalized operating frequency / / f' opt . Plot parameters: /q = c/Xq = 3- 10 8 
Hz, /' = 0.984/o, g/X = 0.05, a/X = 0.08, e r = 60. 
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Figure 7. The: (a) real and (b) the imaginary parts of the moments c p* z and (—m*) as 
functions of the normalized operating frequency / '/ f' opt . Plot parameters: /o = c/Ao = 3- 10 s 
Hz, f' opt = 0.984/q, g/Xa = 0.05, a/A = 0.08. 
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Figure 8. The: (a) real and (b) the imaginary parts of the scattered electric field in the far region 
along the forward direction {cj> = 0), using exact formulas (lim/ C0P ^ +oo E z o tSca t(p, 4> = 0)) 
and the approximate moments model Qim.f. oP ->.+oo E' z0 scat (p, 4> = 0)), as functions of the 
normalized operating frequency / / f' opt . Plot parameters: /o = c/Xq — 3 • 10 s Hz, 
/' = 0.984/o, g/X = 0.05, a/X = 0.08. 



